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Abstract 

In this study, the analytic expressions of the steady probability distribution of tumor cells were 
established based on the steady state solution to the corresponding Fokker-Planck equation. Then, 
the effects of two uncorrelated white noises on tumor cell growth were investigated. It was found 
that the predation rate plays the main role in determining whether or not the noise is favorable 
for tumor growth. 
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I. INTRODUCTION 



As tumors seriously threaten human health, extensive attention has been paid to 
this issue by researchers in various fields [l- 2?]]. It is known that tumor cell growth is 
a complex process, and is governed by environmental fluctuations such as the people's 
spiritual status, as well as other diseases from which they suffer. Recently, researchers from 
the field of nonlinear physics have introduced noise into the model of tumor cell growth, 
where noise refers to the various disturbances involved in tumor growth. For example, Ai 
and coworkers studied the effects of correlated Gaussian white noise in a logistic 

growth model. This model is often as a basic model for cell growth, particularly tumor 
cell growth}^. I29] ]. to describe such growth under ideal conditions without fluctuation. Ai 
and coworkers found that noise during tumor cell growth can induce phase transition, and 
that intensive environmental fluctuations may even cause the extinction of tumor cells. 

□ n 

Furthermore, Zhong and coworkers [3, 126|] investigated the random resonance of tumor 
growth with noise. It was found that the steady distribution probability of tumor growth 
changed from a uni-peak state to a bi-peak state when the intensity of multiplicative noise 
increased. An appropriate intensity of multiplicative noise can destroy the mechanism of 
tumor growth. In contrast, superfluous noise can be beneficial for their growth. Mei and 
coworkers [91] investigated the tumor cell growth model in the presence of correlated noises 
and found that the correlation intensity A and correlation time T play opposite roles in the 
static properties and the state transition of the system. An increase in A can produce a 
smaller mean value of the cell population and slow down the state transition. However, an 
increase in T can produce a larger mean value of the cell population and enhance the state 
transition. 

The abovementioned results demonstrate that tumor growth models with noise are closer 
to the real situation, although most models exhibit some differences from the real process 
of tumor growth. In such studies, researchers have attempted to obtain deeper insights into 
the intrinsic mechanisms of tumor growth and provide new ideas for tumor treatment. 
The noises introduced into the tumor growth model are generally a single noise, correlated 
multiplicative noisey, 5|, etc. It has been shown that different types and numbers of noise 

nnnnnn 

may be operating in the context of tumor growth[2|, |3|, [5|, l8Hld. |26| . and the corresponding 
tumor growth behaviors have also been observed to be different. The random fluctuations 
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introduced in our study are different from those in previous studies. In this paper, the 
effects of two uncorrelated Gaussian white noises on tumor growth will be studied. These 
are the effect of additive noise on the birth rate of tumor cells and that of multiplicative 
noise on the predation rate of anticancer cells. It is shown that the effect of noise on the 
tumor growth is mainly determined by the predation rate of anticancer cells. With changes 
in the parameters, the steady distribution probability of tumor growth changes between a 
single steady state and bi-stable state. 



II. THE DETERMINISTIC MODEL OF TUMOR CELL GROWTH 



Lefever and Garay 30(] studied tumor growth under immune surveillance against cancer 



using the enzyme dynamics model. The model is as follows: 

Normal Cells — ► X, 

X — ► 2X, 

X + E — > E — y E Q + P, 



Here, X, P, E , and E are cancer cells, dead cancer cells, immune cells, and the compounds 
of cancer cells and immune cells, respectively. This model reveals that normal cells can 
transform into cancer cells, and then the cancer cells reproduce, decline, and ultimately die 
out. 

Based on the model by Lefever and Garay, we investigate the Logistic model of Verhulst, 
and only consider the growth of tumor cells and anticancer cells. We assume that tumor 
cells satisfy the following equation[31,32]: 

dX r B x(l-f)-P(X), (1) 



dt \ k B J 

where X is the relative number of tumor cells, is the birth rate of the tumor cells, 
and ks is carrying capacity. P(X) represents predation generated by anticancer cells. We 
take the P(X) expression suggested by Ludwing[31]. In his work P(X) is expressed by 
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BX 2 /{A 2 + X 2 ), where A is a positive constant and B represents the predation rate of the 
anticancer cells. 

As we do not intend to investigate the constant A, it was concealed in the model for the 
convenience of discussion. The transformation parameters can be given by jj], 5, 7|: 

X K t 

x = ~j,r = Ar B ,q = -j-,t = -,(3 = B, (2) 

and by substituting this into Eq. (1), we obtain: 



dx A x\ (3x 2 

— = rx 1 

dr \ q J 1 + x 



where r is the tumor cell growth rate and (3 is the predation rate of anticancer cells. By 
letting f = 0, we can obtain the steady states of the system from the Eq. (3). Clearly, one 
of the solutions is x = 0, while the other solutions satisfy: 

'H)=t£?- (4) 

When r varies and keeps (3 and q constant, the number of solutions (namely equilibria) 
changes between one and three [32]. The range with three solutions changes with the values 
of f3 and q. This also occurs for a variable (or q) and fixed r and q (or (3). 
Based on Eq. (3) (let f(x) = ^), we can draw the curves of f(x) — x as shown in Fig. 
1. For a curve with r = 1.0 and (3 = 2.0, x = and x = x<i are unstable states, since 
df /dx > at x = 0, xi- However, x\ and x% are stable steady states since df /dx < at 
these two points. 

FIG. 1: The equilibria of f(x) vary with a decreases in r and increases in j3 for q = 10.0 (arbitrary 
units). 



In Fig. 1, when r decreases or (3 increases, the solutions 22 and X3 will disappear and only 
the stable state X\ on the left side can be observed (e.g. the curve for r = 1.0, (3 = 3.0). 
In contrast, when r increases or (3 decreases, only the stable state 23 on the right side can 
be observed (as shown in Fig. 2). Clearly, x\ is the refuge equilibrium, while X3 is the 
outbreak equilibrium. From a tumor control point of view, we need to keep the number of 
tumor cells in the refuge state rather than allowing it to reach an outbreak situation. 
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FIG. 2: The equilibria of f(x) vary with an increases in r and decreases in f3 for q = 10.0 (arbitrary 
units). 



When considering the noise, the range of three equilibria is also related to the noises. Simi- 
larly, the equilibrium state changes between 1 and 3 by varying the parameters. Therefore, it 
is necessary to investigate them by using the steady probability distribution function(SPDF). 

III. TUMOR CELL GROWTH MODEL WITH NOISE 

Equation (3) only describes the tumor growth behavior under ideal conditions, without 
fluctuation. When considering a real situation, external environmental disturbances such 
as the individual's state of health, body temperature, and other disease may affect tumor 
growth. In addition, artificial behavior like chemotherapy may also have an effect. Because 
of these external disturbances, the growth rate of tumor cells and the predation rate of 
anticancer cells may vary greatly. In our work, the effects of additive noise on the birth rate 
of tumor cell and multiplicative noise on the predation rate of anticancer cells are further 
considered. Hence the tumor growth equation considering external disturbances can be 
rewritten as: 



where and T(t) are Gaussian white noises. They have the following properties: 



where a and D are the strength of the noises and T(t), respectively. From this, we can 
derive the corresponding Fokker-Planck equation for the evolution of SPDF based on Eq. 
(5) and Eq. (6). The equation is as follows [11]: 





< f (t)f (0 >= 2a5(t - t), < r(t)r(t) >= 2D5(t - t ), 



(6) 



dP(x,t) dA(x)P(x,t) | d 2 B(x)P(x,t) 
dt dx dx 2 



(7) 



where P(x,t) is the probability of the relative numbers of tumor cells, and: 
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B(x) 



x 



1 + x 2 



+ D. 



(9) 



IV. STEADY STATE ANALYSIS OF THE MODEL 

Usually, what we are concerned with is the steady state. For Eq. (7), when 
we can obtain the SPDF of the tumor cells [11]: 



dP(x,t) 

at 



Pst(x) = 



N 
B~{x) 



exp 



x A(x) 



dx 



or 



PJx) 



N 



where N is the normalization constant, 



exp[M(x)], 



B(x) 

When considering Eqs. (8), (9), and (12) together, we can obtain: 

U(x 



M(x) = ax + bx 2 + cx 6 + d\n\E(x)\ + /In 



V(x) 



+ In 



E(x) 



1 + x 2 



+m arctan H(x) + n (arctan K(x) + arctan L(x)) 



where 



a = 



2ra + qp(D + a) 
q(D + a) 2 ' 
r 

' = 2{D + a) ' 

r 

TO 

= 2{D + af 1 



ra (-3D + a + 2yjD(D + a)) + q(5(D + a) (a — D + \j D(D + a 



Aq(D + o)i^D(D + a) - 2D 
r^(a - D) 



m 



2^(D + a) 2 ' 



0, 



(10) 



(11) 



(12) 



(13) 

(14) 
(15) 
(16) 
(17) 

(18) 
(19) 
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ro (3D - a + 2JD(D + a)) + qf3{D + a) (D - a + J D(D + a)) 
n = ' v -, (20) 



2q(D + a)*^2y/D(D + a) + 2D 
E(x) = D + 2Dx 2 + (D + a)x\ (21) 



U(x) = ^D + \ 2JD(D + a) - 2Dx + J (D + a)x 2 , (22) 



V(x) = V15 - \J 2^D(D + a)- 2Dx + \J (D + a)x 2 , (23) 



2 v / STax + \ 2Jd(D + <j) - 2D 
K(x) = V V , (25) 

I2JD(D + a) + 2D 



2VD + ax - \ 2JD(D + a) - 2D 
L{x) = V V . (26) 

l2y/D(D + a) + 2D 



Equation (11) is the main result in this study. Based on this equation, we can plot the 
figures and obtain the curves for P st (x), r, D, (3, and a. In this way, we can determine the 
mechanisms of tumor growth model. 



V. RESULTS AND DISCUSSION 

A. The relationships between r, (3, and x under SPDF extremum condition 

In order to discuss the effects of the fluctuation on the steady probability distribution 
(SPD) of tumor growth, it is necessary to discuss the relationships between r, j3, and x. The 
condition to obtain the extremum of SPDF is: 

A(x) - B'(x) = 0, (27) 

By considering Eqs. (8), (9), and (27) together, we can obtain: 

/ x\ /3x 2x 2 . . 

r { 1 --,)-TT^-"WT^- - (28) 
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Based on Eq. (28), we can draw the curves of r - x and (3 - x in Fig. 3 and 4, respectively. 
As shown in Fig. 3, when 0.1 < a < 1.0, one value of r (0.9 < r < 1.3) corresponds to three 
values of x with a fixed (3 and q. That is, the SPDF of tumor growth will exhibit two peaks 
(corresponding to the three solutions, namely the two stable steady states in the section 2) 
if r is in that range. For different values of o (or (3), the curves almost overlap with each 
other, except in the position around x — 1. This demonstrates that only when x is around 
1 can we observe the difference of r for different values of a (or /?). 

FIG. 3: Plot of the SPDF extrema as a function of x for different a values, using (3 = 2.25 and 
q = 10.0 (arbitrary units). 



FIG. 4: Plot of the predation rate (3 of anticancer cells as a function of x for different values of r 
and a with q = 10.0 under the SPDF extremum condition (arbitrary units). 



Similarly, adopting r = 1.0 in Fig. 4, one value of (3 (1.7 < [3 < 2.5) corresponds to three 
values of x. The range of the two peaks is different with different r values. 
As Eq. (28) is irrelevant to D, the ranges of the two peaks is mainly determined by r and 
(3. In addition to the aforementioned ranges, for smaller or larger values of r (or (3), each 
curve has only one peak (mono-stability). This corresponds to the situation shown in the 
section 2, but there are four parameters here(r, /3, q and a). 

It is clear that the SPD of tumor growth switches between bi-stability and mono-stability 
when the parameter r changes. As shown in Fig. 5, the positions of the peaks vary with 
the different r values. For a small r value, there is only one peak on the left side, which 
represents a small quantity of tumor cells in healthy people or the annihilation of tumor 
cells. With the increase of r value, two peaks can be observed. With a further increase in 
the r value, the peak number is again reduced to one and its position shifts to the right 
side, indicating the steady growth of the tumor cells. Comparing the growth behaviors 
under different r values, it can be concluded that a large value of r is favorable for tumor 
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cell growth. 



FIG. 5: Plot of P s t(x) against x for different r values, using /3 = 2.3, q = 10.0, D = 0.5, and a = 0.5 
(arbitrary units). 



FIG. 6: Plot of P s t(x) against x for different (3 values, using q = 10.0, r = 1.0, D = 0.5, and a = 0.5 
(arbitrary units). 

Fig. 6 shows the effects of the /3 value on the tumor growth. With the increase of the (3 
value, the number of peaks on the curve changes from one to two, and finally back to one. 
It is clear that a large predation rate is unfavorable for tumor growth. 

B. The effects of fluctuations on tumor growth rate 

Figure 7 shows the effects of D on tumor growth for a small ft value (e.g. (3 = 1.7). 
When the value of D increases, the peak intensity of the probability density decreases, 
which demonstrates that an increase in D is unfavorable for tumor growth. 
In the case of a large f3 value (e.g. f3 = 2.6), the peak position shifts to the left side, as shown 

FIG. 7: Plot of P 8 t(x) against x for different D and a values, using (3 = 1.7, q = 10.0, and r = 1.0 
(arbitrary units). 

in Fig. 8. When the D value increases, the SPD of tumor cells moves to the more positive 
x direction. When healthy people receive chemotherapy, normal tissue cells and anticancer 
cells may be killed rather than tumor cells. Therefore, in such a case, chemotherapy is 
favorable for tumor growth (although there is usually a small quantity of tumor cells in 
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healthy people, they can be controlled by the immunity of the human body). 
It can thus be concluded that (3 plays a very important role in tumor growth, especially 
if (3 is very small or very large, while a is less important for the tumor growth. In a real 
situation, the magnitude of the predation rate can determine the body's anticancer ability. 
For a moderate (3 value, for example (3 = 2.26, there are two peaks for all the D values, as 

FIG. 8: Plot of P s t{x) against x for different D and a values, using (3 = 2.6, q = 10.0, and r = 1.0 
(arbitrary units). 

shown in Fig. 9 and 10. For a fixed /3, the SPDF curve becomes flatter with an increase in 
D. The noise interferes with the tumor growth as well as the predation rate of anticancer 
cells. For example, by using chemotherapy, tumor cells may be extinguished. At the same 
time, normal tissue cells and anticancer cells may also be damaged. Here, whether or not 
the increase of D is favorable for tumor growth is determined by a. For a small value of a, 
as shown in Fig. 9 (a — 0.1), the intensity of the right peak decreases with an increase in 
D, which means that it is unfavorable for tumor growth. In contrast, for a large value of 
a, as shown in Fig. 10 (a = 0.8), even though the intensities of two peaks decrease at the 
same time, the SPD of tumor cells moves to the more positive direction of x. That is, when 
a is larger, interference with the anticancer cells is dominant. Thus, the increase of D in 
the case of a larger a when (3 is a middle value is favorable for tumor growth. 
For a different value of r, for example, r = 1.2, the results indicate that the laws acting on 

FIG. 9: Plot of P s t{x) against x for different D values, using (3 = 2.26, q = 10.0, r = 1.0, and 
a = 0.1 (arbitrary units). 

FIG. 10: Plot of P s t(x) against x for different D values, using (3 = 2.26, q = 10.0, r = 1.0, and 
a = 0.8 (arbitrary units). 

Pst(x) and D are same as that of r = 1.0. However, the range of the two peaks is different. 



10 



If the curves for r = 1.0 shift to the more positive direction of (3, they can nearly overlap 
with the curves for r = 1.2. For example, Fig. 11 is almost same to Fig. 9, except that 
the relevant (5 value (2.72) is higher by around 0.46 larger than that for r = 1.0 (ft — 2.26). 
This means that the laws of tumor growth are similar even with a different tumor growth 
rate. 

FIG. 11: Plot of P s t{x) against x for different D values, using f3 = 2.72, q = 10.0, r = 1.2, and 
a = 0.1 (arbitrary units). 



C. The effects of fluctuations on predation rate 

Different values of r, /3, D, and a are adopted in this study. The corresponding results 
show that the effects of fluctuations on j3 are similar to those on r. For a small value 
of j3, an increase in a is unfavorable for tumor growth. For a large value of j3, on the 
other hand, an increase in a is favorable for tumor growth. This effect is independent 
of D. However, for a moderate value of /3, whether or not the noise is favorable for 
tumor growth is determined by D. For a small value of D, interference with (3 is 
dominant; this is favorable for tumor growth. In contrast, large values of D are unfavor- 
able for tumor growth. If r is changed, the laws are same as those mentioned above, but 
the corresponding (5 will be different. This is also similar to the effects of the fluctuation of r. 

VI. SUMMARY 

We investigated the effects of the environmental disturbances on tumor cell growth. By 
solving the corresponding Fokker-Planck equation, we obtained analytic expressions of the 
steady state probability distribution of tumor cells. It was found that the effects of noise 
on the tumor growth are mainly determined by the predation rate (3: (1) For a small value 
of (3, the effects of the disturbance on tumor growth and anticancer cells are unfavorable 
for tumor growth; (2) A large value of f3 is favorable for tumor growth. (3) For a moderate 
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value of P, the effect is determined by the fluctuation in the relative strength of the two 
noises: (a) If the fluctuation strength of the predation rate a is small, the increase of the 
tumor growth rate fluctuation intensity D is unfavorable for tumor growth; in contrast, it 
is favorable for tumor growth if a is large; (b) if D is small, increasing a is favorable for the 
tumor growth; in contrast, decreasing a is unfavorable for tumor growth. Although further 
work is still necessary, it is believed that the present results can give some useful insights 
for the clinical treatment of tumors. 

This work was supported by National Natural Science Foundation of China under Grant 
No. 30600122. 
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